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Abstract 

Quantum dynamics of a general dissipative system investigated by 
its coupling to a Klein-Gordon type field as the environment by in- 
troducing a minimal coupling method. As an example, the quantum 
dynamics of a damped three dimensional harmonic oscillator inves- 
tigated and some transition probabilities indicating the way energy 
flows between the subsystems obtained. The quantum dynamics of a 
dissipative two level system considered. 

1 Introduction 

In classical mechanics dissipation can be taken into account by introducing a 
velocity dependent damping term into the equation of motion. Such an ap- 
proach is no longer possible in quantum mechanics where a time-independent 
Hamiltonian implies energy conservation and accordingly we can not find a 
unitary time evolution operator for both states and observable quantities 
consistently. 

To investigate the quantum mechanical description of dissipating systems, 
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there are some treatments, one can consider the interaction between two sys- 
tems via an irreversible energy flow [1,2], or take a phenomenological treat- 
ment for a time dependent Hamiltonian which describes damped oscillations, 
here we can refer the interested reader to Caldirola-Kanai Hamiltonian for a 
damped harmonic oscillator [3]. 

H(t) = e- 2 ^f- + e^^Wg 2 . (1) 

Zi lit Zi 

There are significant difficulties about the quantum mechanical solutions of 
the Caldirola-Kanai Hamiltonian, for example quantizing such a Hamilto- 
nian violates the uncertainty relations or canonical commutation rules and 
the uncertainty relations vanish as time tends to infinity. [4,5,6,7,8] 
In 1931, Bateman [9] presented the mirror-image Hamiltonian which con- 
sists of two different oscillator, where one of them represents the main one- 
dimensional damped harmonic oscillator. Energy dissipated by the main 
oscillator completely will be absorbed by the other oscillator and thus the 
energy of the total system is conserved. Bateman Hamiltonian is given by 

H = — + Jr—(xp — xp) + (k — j— )xx, (2) 
m 2m 4m 

with the corresponding Lagrangian 

L = mix + — (xx — xx) — kxx, (3) 

canonical momenta for this dual system can be obtained from this Lagrangian 

as 

dL ± p_ dL . (3 

p = — = mx - -x, p = — = mx + -x, (4) 

ox 2 ox 2 

dynamical variables x,p and p,x shoud satisfy the commutation relations 

[x,p] = i, [x,p]=i, (5) 

however the time-dependent uncertainty products obtained in this way, van- 
ishes as time tends to infinity. [10] 

Caldirola [3,11] developed a generalized quantum theory of a linear dissipa- 
tive system in 1941 : equation of motion of a single particle subjected to a 
generalized non conservative force Q can be written as 

d , dT dT 8V , 
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where Q r = —f3(t)J2 a rjQj, and a rj 's are some constants, changing the vari- 
able t to t*, using the following nonlinear transformation 

t* = X (t), dt = cf>(t)dt*, 0(t) = e /oW, (7) 

together with the definitions 

9* = ^, L* = L(q,q*,n P* = |p (8) 



the Lagragian equations, can be obtained from 

d ,dL\ dU 



dt* dq* dq 



0. (9) 



where H* = J2p*Q* ~ L*. Canonical commutation rule and Schrodinger 
equation in this formalism are 

[q,P m ]=i, H ** = i W> (10) 

but unfortunately uncertainty relations vanish as time goes to infinity. [10] 
Perhaps one of the effective approaches in quantum mechanics of dissipa- 
tive systems is the idea of considering an environment coupled to the main 
system and doing calculations for the total system but at last for obtaining 
observables related to the main system, the environment degrees of freedom 
must be eliminated. The interested reader is referred to the Caldeira-Legget 
model [12,13]. In this model the dissipative system is coupled with an envi- 
ronment made by a collection of N harmonic oscillators with masses m n and 
frequencies u n , the interaction term in Hamiltonian is as follows 

N N 2 

tf'—sEc^ + g 2 ^— (11) 

n=l n=l Z " L n^ n 

where q and x n denote coordinates of system and environment respectively 
and the constants Cfi clX6 called coupling constants. 

The above coupling is not suitable for dissipative systems containing a dis- 
sipative term proportional to velocity. In fact with above coupling we can 
not obtain a Heisenberg equation like q + u> 2 q + j3q = for a damped har- 
monic oscillator, consistently. In this paper we generalize the Caldeira-Legget 
model to an environment with continuous degrees of freedom by a coupling 
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similar to the coupling between a charged particle and the electromagnetic 
field known as the minimal coupling. In section 2, the idea of a generalized 
minimal coupling method is introduced and in section 3 the quantum dy- 
namics of a three dimensional damped harmonic oscillator is investigated. 
In section 4, we study quantum dynamics of a dissipative two level system. 
Finally in section 5, dynamics of quantum field of reservoir is investigated. 



2 Quantum dynamics of a dissipative system 

Quantum mechanics of a dissipative system can be investigated by introduc- 
ing a reservoir or an environment that interacts with the system through a 
new kind of minimal coupling term. For this purpose let the damped sys- 
tem be a particle with mass m under an external potential v(x). The total 
Hamiltonian, i.e., system plus environment, is 

H = &~® 2 +v(x)+H B , (12) 

where x and p are position and canonical conjugate momentum operators of 
the particle respectively and satisfy the canonical commutation rule 

[x,p\=i, (13) 

H B is the reservoir Hamiltonian defined by 

/+oo _, 
dPku&lWbiit), u s =\k\. (14) 

-oo 

Annihilation and creation operators b%, bt, in any instant of time, satisfy the 
following commutation relations 

h(t),bt,(t)] = t(k-k>), (15) 

and we will show later that reservoir is a Klein-Gordon type equation with 
a source term. Operator R have the basic role in interaction between the 
system and reservoir and is defined by 

-> r+oo _, 

R(t) = / d^k\f{u^{t) + />£)&£(*)]*, (16) 
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let us call the function f(oJ%), the coupling function. It can be shown easily 
that Heisenberg equation for x(t) and p leads to 

1* -rrr -»l P ~ R 

x = i\H, x\ = , 

m 

p=i[H,p\ = -Vv, (17) 
where after omitting p, gives the following equation for the damped system 

mx*=-Vv-R. (18) 
Using (JTHJ) the Heisenberg equation for bk, is 

k = %] = -tw^g + tr^k.S, (19) 
with the following formal solution 

= %(0)e^ + f dt'e-^^m, (20) 

w 

substituting bAt) from (J2Uj) into (fTHj) . one obtains 
mf + /* dt'x*{t')Mt - t') = -Vv + f (i) 

JO 

gyj- roo 

e(*) = < / ^(/(^ft^e-*** - /*K^ (0)e^)t (21) 

J-oo K 

It is clear that the expectation value of £(t) in any eigenstate of Hb, is zero. 
For the following special choice of coupling function 



\ Att 2 u 



k 



equation (J2TJ) takes the form 

mx + (3x = —Vv + 



i(t) = iJ^L £j jt(b n (0)e-^ - &t( )e*tf). (23) 



In the following we investigate the quantum dynamics of a two level system 
and an isotropic three dimensional harmonic oscillator, both interacting with 
a reservoir, as prototypes of dissipative models. 
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3 Quantum dynamics of a three dimensional 
damped harmonic oscillator 

3.1 Quantum dynamics 

For a three dimensional harmonic oscillator with mass m and frequency u 
we have v(x) = ^rrP^x 2 and therefore we can write (|23|) as 

S+ls + u 2 x = (24) 
m m 

with the following solution 

x(t) = e-^(Ae iuJlt + Be~ luJlt ) + M(t), 

-, r+co 

M(t) = i d 3 k 

J — oo 



JLf_W_ e ^ _ 4(°) ^ 

(25) 



where o>i = yu; 2 — Operators A and £?, are specified by initial condi- 
tions 



A + B = f(0)-M(0) 

f (o) - m(o; 

#0) - R(0) 



:^+^i)I+(^-^)i = f(o)-M(o) 



m 



-M(0), (26) 



solving above equations and substituting A and S in (|25|). one obtains 

af(t) = e~^{ ^-^- sin^xt + x(0) cosa;it H — — x(0) sin^t 

mwi 2ma;i 

fi( ^ sina;it — M(0) sin wit - M(0) cosa;it - — M(0) sina;it} + M(t), 



mu;i 2mwi cji 

(27) 
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also substituting x(t) from (|27)l in (J2U|) we can find a stable solution for b%(t) 
in t —>■ oo as 



3(3 e-**# 2^,. m-m 

-{UJ X{{J) + lUJj, 



\ k V k rn K' 



Hit) = %(0)e-^ - i. 

-M(0)u; 2 -iu^M(O)} 

3/?i r+™ d 3 k' b g {0) sin ^T^' f 



fe V «' k' m 



k' m k' 2 



now substituting for(i) from (J28|) in (JTfiJl and using (|T7jl. one obtains p = 
mx + R. 

A vector in fock space of reservoir is a linear combination of basis vectors 
\N(h) 7 N(k 2 ), ...) B = Y = fc ^ "• |0) B (29) 

where are eigenstates of an d the operators and &t act on them as 



bj;\N(h), N(k 2 ), ...N(k), ...) B = \[N$)\N{k 1 ),N{k 2 ), ...N{k) - 1, ...) B 

bt\N(h), N(k 2 ), ...N(k), ...) B = 0NT® + l|iV(fci), iV(fc 2 ), ...iV(fc) + 1, ...) B 

(30) 

If the state of system in t = is taken to be \i/j(0)) = \0)b <8> |^i, ^2, ^3)^ 
where |0) # is the vacuum state of the reservoir and \ni,n 2 , n 3 ) w is an excited 
state of the Hamiltonian H$ = f — h ^mco 2 ^ 2 , then it is clear that 



(-0(0)1^(0)1-0(0)) = (m + n 2 + n 3 + ^)u. (31) 



7 



On the other hand from (|27|). (|28|) and we find 

/zm^ooK^O)! : imf 2 + ^mu 2 x 2 : |^(0))] = 0, 

^77^4(^(0)1 : H.(t) : |-0(O)>] 
/? 2 w 4 . /"+ 00 da; e ixt 



/im^ool / (ni,n 2 ,7i 3 |f (0)|ni, n 2 , n 3 ) 

27T 2 777 J-oo X CJ 2 - X 2 + l~ X 

m 

f3 2 (3 2 3 

- 7T~ (^l,r72,^3^ 2 (0)|77i,77 2 ,n 3 ) = — — (llj + 77 2 + 77 3 + -), (32) 

2m 2m A u 2 

where : : denotes the normal ordering operator. Now by substituting 
b%(t) from (|2"KJ) into ((H)) , it is easy to show that 



Hm t ^4(#))|:# B (t):h/;(0))] 



, x 77 1 ,77 2 ,r7 3 |f 2 (0)|77 1 ,772,773) 

7T Jo (a; 2 - x 2 ) 2 + ^x 2 

/■^ /* *^ doc i I o / \ i \ 

{ni,n 2 ,n 3 \p {0)\n ll n 2 ,n 3 ) 



nm 2 Jo (uj 2 - x 2 ) 2 + 4j 2 

,3 



/3ar, 3. f°° dx 

(ri! + 77 2 + 77 3 + -) / — ONO ; fl2 - + 



7T777 



2 Jo (uj 2 - x 2 ) 2 + ^jx 



0u . 3. f 00 x 2 dx 

-(77i+r7 2 + 77 3 + -) 



7T777 V 11 Z 2' Jo [to 2 - X 2 ) 2 + ^X 2 

(33) 

For sufficiently weak damping that is when (3 is very small , the integrands 
in (j3l?j) have singularity points x = ±(o;i ± and by using residual calculus 
we find 

^ttw,[(^(0)| : H B : \ifi(0))] = (ni + n 2 + n 3 + ^)co, (34) 

Comparing (|3*T|) and (J53j) . show that the total energy of oscillator has been 
transmited to the reservoir and according to (|32|). the kinetic energy of os- 
cillator tends to zero. 

- H B 

If the state of system in t = is p(0) = Pb ® \S) U where p T B = — - — _u B 

is the Maxwell-Boltzman distribution and \S) U is an arbitrary state of har- 
monic oscillator , then by making use of Ttb [ftt(0)6p {0)pb) = &<y ^ k ^ one can 



eKT -1 
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show that 



lim t -nx>(- -mx (t) H — mu x (t) :) = / m ; dx 

K 2 w 2 W 7 vrm 2 Jo [( w 2 - x 2 ) 2 + 4a; 2 ](e^ - 1) 
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H — -^r / ~m x dx. (35) 

urn? Jo [( a; 2_ a; 2)2 + |i x 2]( e7 ^_ 1 ) 

3.2 Transition probabilities 

We can write the Hamiltonian (|T2"j) as 

H = H + H', 

3 3 
H = Y,( a } a j + 9V + H B> 

i=i 

H' = -^.R+^-, (36) 

where ay, a] j = 1, 2, 3 are annihilation and creation operators of the har- 
monic oscillator. In interaction picture we can write 

ajI (t) = e iHot aj (0)e- iHot = %(0)e- iarf j = 1, 2, 3, 

b Zl {t) = e iHo \(0)e' iHot = %(0)e-^*, (37) 

the terms —.p and ^ are of the first order and second order of damping 
respectively, therefore, for a sufficiently weak damping, is small in com- 
parison with — .p. Furthermore ^- has not any role in those transition prob- 
abilities where initial and final states of harmonic oscillator are different, 
hence we can neglect it in H' . Substituting ajj j = 1,2,3 and b^j from 

(J37J) in — .p, one can obtain Hj in interaction picture as 



e -i(w-wg)t 



-/(^)aA-(0)e- i( ^+-)* + /*( Wfc )at6j(0)e i ( w « + «)*]A; i , (38) 

where fcj j = 1, 2, 3 are the cartesian companent of the vector /c. The terms 
containing just ajb^(0) and aj&t(0) violate the conservation of energy in the 
first order perturbation, because ajb^O) destroys an excited state of harmonic 
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oscillator while at the same time destroying a reservoir excitation state and 
Ojbt(0) creates an excited state of harmonic oscillator, while creating an 
excited reservoir state at the same time, therefore we neglect them due to 
energy conservation and write Hj as 

I ~t i r H~oo 

H> = -i^E j ^ d*k[f{u> % )a%{U)e^-^ - r(^0)e^. 

(39) 

The time evolution of density operator in interaction picture is [14] 

pi (t) = U I (t ) t ) Pl (t ) C/J (Mo) , (40) 
where £// is the time evolution operator, which in first order perturbation is 



UAt, t = 0) = 1 - % [ dtiH'AU) 
Jo 



„• 1 ^ — oo 



2m . =1 . 



2 

Let p/(0) = 1^1,^2)^3)0; w(wi,n2,ri3| <S> |0)b b(0| where \0)b is the vacuum 
state of the reservoir and \nx,n2,ns) u , is an excited state of the harmonic 
oscillator, substituting Uj(t, 0) from (jlT]) in pOjl and tracing out the reservoir 
parameters, we obtain the probability that the reservoir absorbs energy u 
from oscillator 

r (ni+n2+n3+ | W(ni+n2+n3+ i )w = Tr s [{\n x - l,n 2 ,n?,) w w {\n x - l,n 2 ,n 3 \ + 
+\ni,n 2 - l,n 3 ) w u,(|ni,n 2 - l,n 3 | + K,n 2 ,n 3 - l) u w (|rii,n 2 ,% - l|)p s/ 
_ 27rcj(n! + n 2 + n 3 ) f+°° 4 2 sin 2 ^p^ 

3^ y ( ^ )2 > 

where Tr s denotes taking trace over harmonic oscillator and p s i{t) is obtained 
by taking trace of pi(t) over reservoir parameters i,e p s i(t) = Trs^piit)) and 
we have used the formula TtbW^b B^-gW = S(k — k'). 

sin 2 ("p~") f 

For very large times, we can write — ^ — = 2irt5(ujf— u) which leads to 

) 2 



4 7 r 2 u; 5 (n 1 +n 2 + n3)t|/(cu)| 2 
3m 



r _ w V'"l t /6 2 -r y^;i 

1 (ni+n 2 +n3 + |)w^(n 1 +n2+n 3 + i)a; ~ „„ 1 l 4C V 
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For the special choice (J22j) . above transition probability becomes 

_ (ni + n 2 + n 3 )pt 

1 (ni+n 2 +n 3 +|)a;^(ni+n2+n3+i)a; ! 1 44 J 

in this case the oscillator can ont absorb energy from reservoir. 
Now consider the case where the reservoir is an excited state in t = for 
example p 7 (0) = |rai,ra2,n 3 ) w u>(ni,n 2 ,n 3 \ ® b(Vi, -Vrl where 

jl^, ...lp> r )B denotes a state of reservoir that contains r quanta with corre- 
sponding momenta p±, ...p r , then 

rr B[4lVi ; -Vr>B s(Vi>-Vrl%'] = 

rr B [%|l pl ,...l p - r ) B B (V u ...VJ^J=^^-p1)5(fc'-p1), (45) 

and from the long time approximation, we find 

_ 47r 2 o; 5 (ni + n 2 + n 3 )t 2 

1 (ni+n 2 +n3+|)w^(ni+n2+n 3 + i)Lj g m \JK U )\ > 

■p 

(rai+n2+«3+f )^^("-i+"2+ri3 + |)^ 



r 



m 1=1 



(46) 



where pn,Pi2,Pi3 are the cartesian components of vector pi . Specially for the 
choice we have 

r (ni + n 2 + n 3 )ffi 

1 (ni+n2+n 3 + |)w^(ni+n,2+rt3+i)a; ~~ J 

T ' 

1 (rai+n2+n 3 + |)ww^(ni+n2+n3 + |)tij 

= £ 5 Kl - W ) [("1 + + ("2 + l)Pfi + ("3 + 1)P? 3 ] 

(47) 

Another important case is when the reservoir has a Maxwell-Boltzman dis- 
tribution, so let pj(0) = 77,2, n3) w u)(ni, n 2 , n 3 \ ® p^ where 
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-H 



B 



Pb = — e k '-h b ; then using the relations 



Tr B [b^p T B bg}=Tr B [blp T B bi} = 0, 



6 KT — 1 



KT 



TrrftfM = 6{k - , (48) 

eifT — 1 

we can obtain the following transition probabilities in very long time 

_ An 2 u 5 (n 1 + n 2 + n 3 )t\ f(u)\ 2 e^r 

_ 47r 2 ^ 5 (m + n 2 + n 3 + 3)t |/(u;)| 2 



3m exr — 1 ' 



(49) 

substituting (J22j) in these recent relations we find 

(ni + n 2 + n 3 )f3te^r 

1 (n 1 +n 2 +n i + l)u)-^{n 1 +n 2 +n s +\)u) ~ m ( e WT — ~T) ' 

(ni + n 2 + r2, 3 + 3)/3t 



■n 

1 (ni+ra 2 +n 3 + §)— (rn+na+ns+f) mie^T - 1 



(50) 



So in very low temperatures the energy flows from oscillator to the reservoir 
by the rate r (ni+n2+n3+ | )a ,^ (ni+n2+fl3+ i )w = feilglZ*^ an d no energy flows 
from the reservoir to the oscillator. 



4 A dissipative two level quantum system 

Let us write the Hamiltonian (|12|) as 

H = H s + H B --.p+^-, (51) 
m 2m 



12 



where H$ = + v(x) is the Hamiltonian of the system. The H s has two 
eigenvalues E\, E 2 corresponding to eigenkets |1), |2) respectively, and we can 
write 

H s = l(E 2 -E 1 )a z + ^{E 1 +E 2 ) a z = |2)(2| - |1)(1|. (52) 

For sufficiently weak damping we can neglect the term because it is 
proportional to the second order of damping and write H as 

H = ~u a z + H B - J d 3 k[(G 12 a + G 21 a% - J d*kb\[{G* 2l a + G* l2 a^] 
Gij = ^^-k.pij = iuJijf(u^)k.Xij 

Xij = (i\x\j) Pij = (i\P\j) i,j = 1,2, (53) 

where u>q = E 2 — E\ and o = |1)(2|. Using commutation relations 

[a,J] = -a z [a,a z ]=2a [a z ,^} = 2a\ (54) 
one can easily obtain Heisenberg equations for the two level system 

cr = i[H, a] = —iuj a — i J d 3 k[G 2 ia z b^ + G{ 2 bla z ] 

& z = i[H,a z ] = -2i J d 3 k[G l2 abj: - G 2l a^b n ] - 2% J d 3 k[G* 21 b\a - G* l2 bta% 

(55) 

Since equal-time system and reservoir operators commute, we can write the 
Heisenberg equations (J55)l in different but equivalent ways. For example 
we can use the normal ordering were annihilation operator of the reservoir 
bj: appears at the right and creation operator bt appears at the left of the 
system operators, i.e, a, o\a z , or we can use antinormal ordering wereh 
annihilation operator of the reservoir b^ appears at the left and creation 
operator &t appears at the right of the system operators, i.e, a,a\a z . In 
equations (|55j) we have used the normal ordering. 
One can easily obtain the Heisenberg equation for fr^ as 

bj- = i[H, y = -iw^ + % [ d 3 kG* 21 a + % f d 3 kG* 12 a\ (56) 
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with the following formal solution 



b i (t)=bz(0)e- iu i t + iGZ 1 / dfa^e'^-^ + iG* 2 dt'aUt^e-^-^. 

Jo Jo 

(57) 

Let us assume that damping is sufficiently weak that we can take the ite- 
grands (|57j) we can take 

(7(0 ^ a(t)e-^ o( *'- t) a\t') ^ ^(t)e luJo{t '- t] . (58) 

This is called the Markovian approximation [15] which replaces the system 
operators in (|5T|) by an operator that depends on the system operators at 
the same time t, without taking into account the memory of these operators 
at earlier times. Substituting (J55j) in (|57j) and using (|B3j) we find 



(a) = -iu (a) + J d 3 k\G 21 \ 2 [(a z a) J cft'e^— o)(*'-t) _ ^ J dt > e ^ k+ u> )(t> -t)^ 

+ / d 3 kG 21 G*[(a z a^) f dt'e^^'^) - ( a ^a z ) f dt'e^-^''^] 
J Jo Jo 

(& z ) = 2{aa ] ) J d 3 k\G 12 \ 2 [J* dt'e^^'-^ + j^V^^) (*'"')] 

-2(a ] a) [ ^|G 21 | 2 [/ < ^e i( ^)C'- i ) + /"Ve"^-'*)^-*)], 
J Jo Jo 

(59) 

where we have taken expectation values in state \ip(0)) = |0)b<E> IS 1 ) which is 
a tensor product of the vacuum state of the reservoir \0)b and an arbitrray 
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state of the system, \S). In long time approximation, we can write 



d 3 k\G 21 \ 2 / dt'e^-^'-V 
Jo 

Aituj 2 \x 21 \ 2 f°° 



jo 3 Jo + 



o 



(%|/(^)| — + TT^Wg - Wo)] 

o jo ujj: — loq 

d*kG 2X G\ 2 f^X^^/^l^l 2 fcJ^'H 
Jo 3 Jo Uj + ujq 

(60) 

substitution of ( [BUJ) in (JoUj) gives 

(a z ) = -A^a) = -2/x(l + (<r z )) V = 3 ~ 

(a) = -iu (a) + (iAi - iA 2 - fi)(<r) + (zA 2 - iAi - ^((x 1 ) 

. 47r 2 a;g|x 12 | 2 r . I/MIM » 47r 2 w 6 |f 12 | 2 r |/(^)| 2 ^ 

Ai — / diOt A 2 — / dLOr ■ . 

3 Jo wj-Wo 3 Jo u^ + Uq 

(61) 

The solution of the first equation in (jBTj) is 

(a,(t)) = -l + (l + (a 2 (0)))e- 2 ^ (62) 

If (<r 2 (0)) = 1 , i.e , the two level system is initially in the upper state |2), 
it decays to the lower state (a z ) = — 1 with the rate 2/i, and if the system is 
initially in the lower state, it remains for all times in that state. 
By defining F = a + and E = — a, the second equation in (JBTJl can be 
written as 

(F) = iT(E) — 2fi(F) r = w -2A 2 -2A lj 

(E)=iu (F), (63) 
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with the following solutions 



(f) = c ie iQ+t + c 2 e n -\ n ± = i/i± iJtf + uj y 



If (j 2 + ljqT < 0, we conclude that (a) decays with the rate fi. 



(E) — ^-Cie iU+t + pp-C 2 e in "*. (64) 

i Z_l_ \l — 



5 Quantum field of the reservoir 

Let us define the operators Y(x, t) and LTy (x, t) as follows 

/+oo A 3 h _ _ 

y ; = (%(^ + 6t(t)e- ife - 5 ), 
-oo ^2(2^)3^ 

/+oo / 7~jZ , - -. 

oo ^^^^ - ^ } ' (65) 

then using commutation relations (jl5j) . one can show that Y(x, t) and LTy (x, t), 
satisfy the equal time commutation relations 

\Y(x,t),IL Y $,t)]=i8(x-iP), (66) 
furthermore by substituting b^(t) from (f2T)J) in (|65jl . we obtain 

fE^M = v 2 F + 2f (t).M(f), M(f) = tfe d 3 k^j ^yj(u % )ke 

BY ■ -> -> r+°° f(uA 

U Y (x, t) = — - 2x{t).N{x), N(x) = Im / d 3 k , JK h) =ke~ ik - x , 



Of 



so Y(x, t) satisfies the following source included Klein-Gordon equation 

d 2 Y ■■ ■ 

V 2 Y = 2x(t).N(x) + 2x(t).M(x), (68) 



dt 2 

with the corresponding Lagrangian density 

1 BY 1 -» • <9F • - 

£ = -(— ) 2 - -Vy.VF - 2S.N(x)— + 2x.M(x)Y. (69) 
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Therefore the reservoir is a massless Klein-Gordon field with source 2x.N(x) + 
2x.M(x). The Hamiltonian density for (J58j) is as follows 

K = '- + -|VY| 2 - 2x.MY, (70) 

and equations (joT|) are Heisenberg equations for Y and Ily. If we obtain 
and &t from (J55)l in terms of Y and Hy and substitute them in Hb defined 
in ([Tlf . we find 

r+oo , TT 2 1 -, 

Hb = rf 3 K4% = ~f + ^Y\ 2 - (71) 

6 Concluding remarks 

By generalizing Caldeira-Legget model to a reservoir with continuous de- 
grees of freedom, for example a Klein-Gordon field, a new minimal coupling 
method introduced which can be extended and applied to a large class of dis- 
sipative quantum systems consistently. This method applied to an isotropic 
three dimensional quantum damped harmonic oscillator and a dissipative two 
level system as prototypes of important dissipative models. Some transition 
probabilities explaining the way energy flows between subsystems obtained. 
Choosing different coupling functions in (fTfij) we could investigate another 
classes of dissipative systems. 
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